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An advanced vehicle concept, the HWB N2A-EXTE aircraft design, was tested in NASA 
Langley’s 14- by 22-Foot Subsonic Wind Tunnel to study its acoustic characteristics for var- 
ious propulsion system installation and airframe configurations. A significant upgrade to 
existing data processing systems was implemented, with a focus on portability and a re- 
duction in turnaround time. These requirements were met by updating codes originally 
written for a cluster environment and transferring them to a local workstation while en- 
abling GPU computing. Post-test, additional processing of the time series was required to 
remove transient hydrodynamic gusts from some of the microphone time series. A novel 
automated procedure was developed to analyze and reject contaminated blocks of data, 
under the assumption that the desired acoustic signal of interest was a band-limited sta- 
tionary random process, and of lower variance than the hydrodynamic contamination. The 
procedure is shown to successfully identify and remove contaminated blocks of data and 
retain the desired acoustic signal. Additional corrections to the data, mainly background 
subtraction, shear layer refraction calculations, atmospheric attenuation and microphone 
directivity corrections, were all necessary for initial analysis and noise assessments. These 
were implemented for the post-processing of spectral data, and are shown to behave as 
expected. 


Nomenclature 


AoA 

= Angle of Attack 

B 

= normalized signal bandwidth 

BENS 

= Broadband Engine Noise Simulators 

CSM 

= cross-spectral matrix 

D 

= baffled microphone directivity pattern 

DAS 

= data acquisition system 

EPNL 

= effective perceived noise level 

FRF 

= frequency response function 

G 

= one-sided power spectral density 

HWB 

= Hybrid Wing Body 

J 0 

= Bessel function of the first kind, order () 

K 

= Ivullback-Leibler divergence 

k 

= acoustic wavenumber 

M 

= Mach number 
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N = number of samples in a block of data 

P = probability distribution function 

p = probability density function 

Q = probability distribution function, estimate of P 

q = probability density function, estimate of p 

R = microphone diaphragm radius 

rms = root-nrean-square 

y n = individual sample in a block of data 

a = gamma distribution shape parameter 

/3 = gamma distribution scale parameter 

r = gamma function 

7 = incomplete gamma function 

9 = wave angle of incidence, referenced to microphone diaphragm normal 

v = effective degrees of freedom of a signal of block size N 

a 2 = variance of a block of data 

= sum of the squares of the samples in a block of data 


I. Introduction 

N ASA’s Hybrid Wing Body (HWB) N2A-EXTE acoustic test in Langley Research Center’s 14- by 22- 
Foot Subsonic Wind Tunnel was a large-scale experiment to demonstrate potential technologies to reach 
the Environmentally Responsible Aviation Project’s N+2 system level aircraft noise goal of 42 cumulative 
EPNdB below Stage 4P The test required significant effort in facility modification, facility calibration ] 3 and 
instrumentation.® The outcome of the test supports analysis of turbomachinerj® and jet noisd^ shielding, 
airframe noise and system noise . 7 

Analysis requirements drove the test to involve approximately 2,000 unique arrangements of the model, 
instrumentation and operational parameters (an example is shown in Fig. [Tj) , leading to the acquisition of a 
large amount of aeroacoustic dataP Initial results were desired in near-time, allowing for potential decisions 
regarding additional test configurations. Post-test processing required the development and use of multiple 
types of acoustic corrections to the data to reconstruct best estimates of lossless, free-field measurements. 
These corrections range in complexity from background noise subtraction to analysis and rejection of transient 
hydrodynamic contamination of individual microphone time series. 

This paper presents the data flow of the acoustic processing for the initial analysis of the test results. 
The transient rejection technique is shown to qualitatively reject obvious contamination, while retaining data 
of interest. Individual microphone, spectral-level corrections which are common to most of the initial data 
processing are discussed, with some characteristic results provided. 


II. Initial Data Acquisition &; Processing 


Details of the data acquisition hardware are provided elsewhere EEI In brief, voltage signals from the 
microphones and surface pressure sensors were transmitted through various hardware-appropriate amplifiers, 
filters and power supplies, some controlled automatically by the data acquisition system (DAS) and some 
set manually. These voltage signals were acquired by the DAS, a 192-channel 3-unit assembly of National 
Instruments (NI) PXI-1045 chassis populated with 16-bit NI PXI-6120 digitizers. For any given set of 
model-traverse configurations, the individual channel gains were set such that the observable dynamic range 
was maximized within the limits defined by the 16-bit system, while avoiding over-ranging. These settings 
were matched between data acquisition from noise sources of interest and corresponding background noise 
acquisitions, so any background noise subtraction would maintain proper scaling for both pressure and 


voltage signals (see Section IV.Al. To maximize system throughput, time series data were streamed to a 
RAID system as raw 16-bit integer data from the DAS, with appropriate metadata for conversion from 
integers to volts stored in the header of each time series binary file. 

In-situ data processing had two somewhat competing requirements: accessibility and efficiency. Time 
series stored in a custom binary format allow for efficient disk and memory utilization, but require custom 
code with knowledge of the format, often impeding accessibility. Data stored in, for example, ASCII format, 
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Figure 1. Example arrangement of the HWB model, phased array and tower traverses. 


is generally highly accessible but very expensive in terms of Hie size. Prior to the test, it was decided that 
a common data format tailored for large-scale scientific data was an appropriate archival tool. Specifically, 
the Network Common Data Form (NetCDF) was selectedP A post-acquisition code was written to convert 
the custom binary data to the more accessible NetCDF format. This code was written in a combination of 
National Instruments Lab VIEW®- 10 and C. It stored appropriate DAS information, along with facility and 
model data, in the header and as variables of NetCDF files for every time series, along with the 16-bit integer 
time series data itself. 

In parallel, a code to compute the cross-spectral matrix (CSM) of each acquisition was required. Histori- 
cally this task has been offloaded to a Fortran code operating on a cluster, and has not had a sufficiently-short 
turnaround time to allow for near-time data analysis. The cluster-based Fortran code was replaced with 
a MATLAB implementation operating on a local workstation. The effort of computing the Fourier trans- 
forms and vector outer products of each block of data to construct an average CSM was offloaded to an 
NVIDIA Tesla C1060 GPU installed in the workstation, as GPU computing has native language constructs 
in MATLAB.^ 1 CSMs were stored in NetCDF files with facility, model, and processing information. 

Work Flow Architecture 

The processing workflow in the construction of CSMs from the acquired time series data was as follows 
(illustrated in Fig. [2|: 

1. Load the 16-bit integer time series data. 

2. Convert the data to single-precision floating-point format. 

3. Use the DAS range information to scale the data from integer counts to volts, as observed by the DAS. 

4. Use the amplifier gain settings to scale the data from volts as observed by the DAS to volts output 
from the transducer. 
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5. Use the transducer sensitivity database to scale the time series data from volts to pascals. 

6. Check the transducer type. If it is a microphone flush-mounted in the array plate, divide the pressure 
signal by two to account for pressure doubling on the plate face. 

7. Split the data into blocks of 8192 points. Subtract the block mean, apply a Hamming window, and 
Fourier transform the data. 

8. Apply the appropriate broadband level correction for the window function and correct for the trans- 
ducer frequency response function (FRF). 

9. Compute one-sided auto- and cross-spectra through the traditional RMS averaging procedure!- 12 ! 

Normally, data were acquired at a sampling rate of 250 kHz for just over 30 seconds, yielding 920 blocks of 
data. For a block length of 8192 points, the resultant CSMs had bin widths of 30.5 Hz. 

Implementation 

In practice, many of these steps were combined to maximize throughput. The executed workflow was: 

1. Pre-compute all of the scalar corrections required to convert integer time series data to pressure time 
series data, for every transducer. 

2. Multiply these scalar corrections by a Hamming window for every transducer, and transfer this array 
of scaled window functions to the GPU. 

3. Initialize an all-zeros CSM of complex single-precision on the GPU. 

4. For every block of data: 

• Transfer the 16-bit integer time series blocks for every transducer for the given block number to 
the GPU. 

• Re-cast the data to single- precision floating-point format. 

• Element-by-element multiply the array of blocks by the array of window functions. 

• Batch-FFT the array of data blocks. 

• Compute the outer product of the vector of transducer Fourier coefficients for each frequency and 
add to the initialized CSM. 

5. Transfer the CSM back to workstation memory. 

6. Scale the CSM to account for the window function broadband correction, number of blocks, and 
transducer FRFs. 

The modified process reduced CSM computation time to well under a minute, with most of the time spent 
transferring the data across the control room network. The majority of this improvement was due to efficient 
use of the strengths of GPU computing, which alone accounted for an order-of-magnitude reduction in the 
processing time when compared to an equivalent CPU-based process executed on the workstation. 


III. Post-Test CSM Reconstruction: Transient Analysis 

The in-situ CSMs were treated as approximates. A refined analysis was discussed and implemented post- 
test. Simple changes included an upgrade of the processing workstation with an NVIDIA Kepler K20c GPU, 
implementation of a correction for the Precision Filter FRFs (incorporated at the same stage of processing 
as the transducer FRFs), and an option to compute autospectra only, in order to minimize memory and disk 
storage when cross-spectra are unnecessary. 

A more complex change involved the assessment of the data for transient contamination. NASA Langley’s 
14- by 22-Foot Subsonic Wind Tunnel is, by design, an aerodynamic wind tunnel which can operate in an 
open test section configuration. While significant acoustic improvements have been applied to the facility, 
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Figure 2. Work flow architecture to compute the CSM for a given acquisition. 
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measurement microphones are, under some installation configurations, close enough to the open-jet shear 
layer that hydrodynamic gusts may contaminate the out-of-flow acoustic measurements. This was primarily 
observed when microphones were at the far-downstream end of the test section, although occasional gust 
impingement was seen at other measurement stations. 

An extreme example of this from the airframe noise component of the test is shown in Fig. [3] The plotted 
data are for an acquisition where one of the speakers embedded in the model bodj 2 is driven with a random 
bandpass signal spanning 4 kHz to 16 kHz. The HWB model is pitched to an angle of attack (AoA) of 14.5°, 
and the test section Mach number is M = 0.23. The acoustic hardware is traversed to the far-downstream 
end of the test section. As shown by the time series in Fig. |3^a), the array center microphone signal appears 
as might be expected for a stationary, band-limited random signal. The south tower microphone, located in 
the upper-right-hand corner of the picture in Fig. |T| clearly experiences extreme random bursts as shown 
in Fig. |3jb) . The corresponding autospectra are shown in Figs. |3^c) andj^d). While the array center 
microphone spectrum shows the low frequency content of the signal at 4 kHz, the south tower microphone 
spectrum is masked by the low frequency content of the burst. Note that at this stage of processing, two 
clean signals would not overlay due to differences in propagation distance between the source and each 
microphone, speaker directivity, and additional effects to be discussed later. 

The south tower microphone time series is partitioned into 920 blocks, each of length 8192 points, and a 
histogram of the block variances is computed. This histogram is shown with respect to the left axis in Fig. 
[4j with the 16 most energetic blocks removed from the plot. Even without these blocks, which would extend 
the x-axis beyond a variance of 500 Pa 2 , this histogram shows a long, thin tail in the direction of large 
variance values. While the contamination from high-variance blocks is obvious through visual inspection of 
the time series and can be dealt with by manually evaluating the data and rejecting blocks, applying such 
a procedure to roughly 250,000 time series records is intractable. A novel, automated process for analyzing 
the data is required. 

Several methods were considered for separating the data. Initially, visual inspection suggested that a 
voltage threshold could be used in separating the data. However this was somewhat arbitrary and specific 
to each run condition. A method evaluating histogram shape was then evaluated. Here, data histograms 
of block rms levels were constructed using sub-blocks l/8th the size of the full block length (1024 points), 
allowing for more converged behavior of the histograms. Histograms from known clean microphone signals 
for a given test condition, for example the array center microphone positioned upstream in the test section, 
were also computed. The ratio between the peak sub-block rms and the rms at the histogram peak, or 
the mode rms, was computed to give a sense of the span of the high-variance tail of the histogram. The 
ratio for the reference histogram was compared to the ratio for the contaminated histogram, and sub-blocks 
with rms ratios beyond this reference ratio were flagged. Any full block containing a flagged sub-block 
was rejected. This technique was successful in identifying and eliminating block contamination, although 
the requirement of a reference signal required the assumption that for every run condition an appropriate 
reference was known. A process which avoided the use of a reference signal was considered and selected for 
the final analysis, which follows. 

The core of the selected rejection process is the assumption that the underlying acoustic signal of interest 
is a stationary, zero-mean, gaussian random process. If the acoustic samples, y, are truly gaussian-distributed 
random variables with zero mean and unit variance, then the sum of the squares 

N 

Xn = J2vI, (!) 

n = 1 

of a set of N samples is a random variable itself, following a chi-square distribution with N degrees of 
freedom.^ If both sides of this relation are divided by N, this shows that for N samples, or for data blocks 
of length N, the block variance can be related to the sum of the squares, 



n—1 


( 2 ) 


This also follows a chi-square distribution, by simply normalizing the data values y 2 by N. Note that for 
finite N this relation involved a biased estimator of the block variance. An unbiased estimator can be used 
by substituting N — 1 for N in the normalization. 
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(a) Array center microphone time series. 
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(b) South tower microphone 7 time series. 
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Figure 3. Example data contamination by hydrodynamic impingement. The two compared micro- 
phones observe a calibration signal with an output band of 4 kHz to 16 kHz, emitted by one of the 
model embedded speakers. The HWB model is at 14.5° AoA, and the test section Mach number is 
M = 0.23. Acoustic hardware are at the far downstream traverse location. 


In practice, the acoustic signal is not truly random, but has a finite bandwidth. This normalized band- 
width, B 1 alters the effective degrees of freedom, v, of the signal?-® For example, a block of 8192 samples of 
a signal which is truly random has a spectrum of white noise and a bandwidth of 100%, so u = N = 8192. 
If the signal passes through an ideal lowpass filter set to 50% of the sampling bandwidth, or B = 0.5, the 
effective degrees of freedom is v = B x N = 4096. The fractional, normalized bandwidth can be estimated 
through a simple procedure. First, the one-sided power spectral density of the signal must be computed. 
This function of frequency, G yy (/), must then be normalized such that its peak is unity, 

Gyy„ (/) = max ^ } (/)] ■ ( 3 ) 

The average of this normalized spectral density is then computed by integrating across the measurement 
bandwidth and dividing by the span, 


j dmax 

B = / Gyy^ oim (/) d /. (4) 

/ max J 0 

While the acoustic data sample distribution does have a mean of approximately zero due to the use of a 
high pass filter during data acquisition, the variance of the distribution for y is unknown. A more general 
distribution is necessary to model the block variance, a 2 . As a generalization of the chi-square distribution, 
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Figure 4. Histogram of block variances from the south tower time series in Fig. [ 3 ] excluding the 16 
most energetic blocks, and modeled data probability density functions. 


the gamma distribution can be used !- 15 The probability density function for a gamma distribution of the 
block variance (with a zero location parameter) is given by 


p[a 


l ) = 


— e^ 


/?r (a) V P 

where a = v / 2 is the shape parameter, /3 is the scale parameter and T is the gamma function 


/» OO 

T (a) = / f a_ 1 e _t df. 
Jo 


( 5 ) 


(6) 


For /3 = 2 this collapses to the chi-square distribution. This scale parameter allows a distribution fit to 
handle non-unity variance of y. 

As previously stated, the effective degrees of freedom and thus the shape parameter of a distribution fit 
can be derived from the signal bandwidth, which can be approximated from a spectral estimate of the signal. 
The scale parameter must now be determined. An easy, if biasedf 16 ^ estimate of fi can be obtained from its 
maximum likelihood estimator 

t2 

( 7 ) 


a 


where a 2 is an estimate of the mean of the block variances. When evaluating the mean of the block variances, 
however, it is readily observable that this mean will be sensitive to the high- variance, hydrodynamically- 
contaminated blocks. A statistic less sensitive to extreme values is necessary for computing /?. One such 
statistic is the median of the block variances. The median occurs where the probability distribution function 
is 0.5. The probability distribution function for the gamma distribution is given by 


Pa 2 ) = 


7 (a, i ) 

r» ! 

where 7 is the (non-normalized) incomplete gamma function 15 


( 8 ) 


7 a, 


The equation for the median variance is thus 


t a ~ 1 e~ t dt. 


1 7 


0 ^) 


T(a) 


( 9 ) 


(10) 
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Software libraries exist for efficiently inverting 7 for a given a, thus yielding an estimate of the median 
variance normalized by fi. The experimental median variance can then be divided by this estimate, yielding 
an estimate of fi. Thus for a given shape factor a, two scale factors can be readily computed from the data. 
One, /3 m ean ; is based on the mean of the block variances and susceptible to extreme values of contamination. 
The other, /3 me dian, is based on the median of the block variances, and is less susceptible to extreme values. 

Having two scale factors allows for the construction of two gamma distributions. These can be compared 
to gain some sense of the relative influence of extreme block variances on the data set. Visually this is done in 
Fig. [4]with respect to the right axis. One numeric metric for comparing distributions is the Kullback-Leibler 
divergence.^ The divergence if is a measure of the information lost when probability distribution Q (or 
density q) is used to estimate distribution P (or density p ) . This is expressed as 


K ip Ik) 




p(<7 2 ) d (<T 2 ) . 


( 11 ) 


While in general this can be difficult to compute, it is greatly simplified in the case of two gamma distributions 
with a common a. In this case, some manipulation yields 


K (Pll9) a -q -0. =Q ln /3p + 




or as used in this application 


A iPmedian\\Pmean) c 


„=a = a In A 


In fimedia 


fi median fin 
fimean 


( 12 ) 


(13) 


To summarize, two data distributions can be estimated. The distribution based on the block variance 
mean is more sensitive to high-variance contamination than the distribution based on block variance median. 
A metric is constructed for comparing the two distributions. Now a procedure is proposed for determining 
which blocks of a given time series to retain, and which to reject. It should be noted here that for the number 
of blocks used, converged data distributions are not expected. The intent of the following procedure is to 
provide an automated engineering tool to locate contaminated blocks in the time series, not to accurately 
estimate the probability distribution of the acoustic data block variance. 

First, every microphone time series is broken into 920 blocks of desired length N = 8192 points. The 
variance of each block is computed, and the blocks are sorted by their variance, from low to high. The initial, 
100 lowest-variance blocks are selected. An autospectral density is computed from this set of low-variance 
blocks, and used to compute a. /? mea n and /3 me dian are computed, followed by K. The next block, in order 
of ascending variance, is added to the set of blocks and the process repeated. This continues until all of the 
blocks of data have been used, producing K as a function of block count in order of ascending variance. The 
block count of minimum \K\ is selected as the cut-off. Sorted blocks above this are rejected. 

As stated, 100 blocks comprise the initial lowest-variance block set. This minimum provides a safeguard 
should the automated algorithm attempt to reject too many blocks. When computing isolated autospectra, 
each time series is handled independently with the 100 block minimum. For CSM construction, if a block is 
rejected for one channel of the system it is rejected for all channels of the system. This may allow a CSM 
to drop below 300 blocks, which based on prior experience has been suggested as a minimum requirement 
for converged cross-spectral data. If this is the case, the CSM is considered sufficiently contaminated that 
the rejection procedure may not improve the results. For such a situation (highly infrequent), the rejection 
technique is disabled, the condition logged in the CSM file header, and all data blocks are used. 

The transient rejection procedure is applied to the south tower microphone data from Fig. [|;b). Of the 
920 input blocks, 567 are rejected. The histogram of the remaining block variances is shown in Fig. [5j along 
with the median- and mean-based probability density functions for the retained block set. As expected, the 
probability density functions overlay for the minimum value of \K\. The output of the procedure is shown in 
Figs. [6](a) and[6](b). Visually, the technique has identified and removed the obvious contamination from the 
time series. In the spectral analysis, the 4 kHz content of the signal is now visible, with a reduction of up to 
10 dB in the microphone autospectrum at lower frequencies. Higher frequencies are unaffected. These results 
indicate that transient signal analysis and rejection may be critical in the use of potentially-contaminated 
data sets. 
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Figure 5. Post-rejection histogram of south tower time series data from Fig. [ 4 J along with post- 
rejection models (mean model near-completely overlays the median model). 
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(a) Blocks accepted and rejected by the algorithm. 



Figure 6. Results of transient rejection algorithm when applied to the south tower time series data 
from Fig. [3](b). Data blocks are plotted as a function of time. The shift in the estimated data 
autospectrum is shown. 
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IV. Post-Test Data Corrections 


A series of data corrections are often necessary prior to application and interpretation of the autospectral 
data from the HWB wind tunnel test. All of these corrections are applied on a narrowband basis. In 
autospectral analysis, all of the corrections aside from background subtraction assume a well-defined acoustic 
source location. At frequencies where this is not the case, either due to ill-defined acoustic source regions or 
the dominance of extraneous noise, these techniques would be unreliable. 

For beamforming and deconvolution analysis the procedure is similar. However, every correction which 
must assume a source location instead generates a set of corrections for every combination of beam map grid 
points and microphones, and for every frequency bin. 

IV. A. Background Subtraction 

Background subtraction is intended to remove extraneous noise from that of the signal of interest. It is 
performed without regard to an acoustic source location. Background subtraction is a simple data correction, 
which has been used repeatedly for previous NASA Langley aeroacoustic experiments.®^® Here, auto- and 
cross-spectra of a facility measurement without the acoustic source of interest are subtracted vectorially from 
auto- and cross-spectra of a facility measurement with the acoustic source of interest. As stated previously, 
this is done for runs with matched acquisition settings (identical ranges and gains), so both voltage and 
pressure noise are scaled properly between the two measurements. Subtraction occurs in the power domain 
with units Pa 2 for both the real and, for the CSM off-diagonal cross-spectra, imaginary components of the 
spectra. 

In application, normally background subtraction assumes that there is no correlation between the back- 
ground noise measurement and the source of interest, so any strong interaction effects may interfere with this 
procedure. This correction is the first of the post-processing corrections applied, as it makes no assumption 
about the nature of the acoustic field beyond the independence of the background noise from the source of 
interest. Often in low signal-to-noise ratio frequency ranges, the background noise level for certain frequency 
bins for the autospectra may be equal to or greater than the noise from the source of interest. In this case, 
the power for the bin is often set to zero in applications, the bin flagged for user notification, and the bin not 
considered in subsequent analysis. At user discretion for spectra where this happens to a significant number 
of bins in frequency bands of interest, background subtraction may be excluded from the data processing. For 
the cross-spectral terms, there is no distinct error determination measure, such as positive or negative sign. 
Applications determine that. For example, CSM-based array beamforming may be substantially immune 
from low signal-to-noise issues. 

A simple example of background subtraction for autospectra is illustrated in Fig. [TJ Here, the acoustic 
hardware were traversed to the far-downstream location and the tunnel operated at M = 0.17. Model 
airframe noise, along with the tunnel wind noise, were included in the data acquired. A model-embedded 
speaker was then turned on, and data were acquired again. Background subtraction extends the range of 
observable speaker input signal at lower frequencies. At higher frequencies, it extends the observable range of 
the speaker output. Note that the speaker input range was 4 to 16 kHz, but system non-linearities extended 
the noise output range. Power bins which are set to zero, near 3 kHz, are not plotted. Near-zero powers 
such as the trough near 37 kHz are seen. 

IV. B. Shear Layer Correction 

In open-jet test section aeroacoustic wind tunnel testing, shear layer correction is necessary for analysis of 
the data. While the overall level correction to the data may be minor at low Mach numbers, the shear layer 
plays a significant role in the path traveled by a packet of acoustic energy. For this study, Amiet’s shear 
layer correction method is used. 21 Here, the facility shear layer sides are modeled as planes, and an acoustic 
source location must be defined. This source location is determined by the particular test configuration. The 
shear layer correction technique determines the transmission path for a packet of acoustic energy traveling 
from the source to an observer. Additionally it calculates the emission coordinates for the source-observer 
position, which follow the phase fronts of the acoustic waves emitted from the source and correspond to 
the emission angles from an in-flight source passing an observer. An example of a shear layer analysis for 
a test case dominated by the main landing gear airframe noise is shown in Fig. [8] This figure shows the 
correction for a microphone located above the model. From Fig. |T] it is evident that there are three planar 
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Figure 7. Example of background subtraction process applied to the array center microphone observ- 
ing a model embedded speaker. The speaker signal has an output band of 4 kHz to 16 kHz. The 
HWB model is at 14.5° AoA, and the test section Mach number is M = 0.17. Acoustic hardware are 
at the far downstream traverse location. 



Figure 8. Example of shear layer correction applied to a microphone stationed 365 inches downstream 
from the test section entrance. The HWB model is at 14.5° AoA, and the test section Mach number 
is M = 0.17. 
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free shear layers bounding the test section. The analysis procedure must determine which shear layer refracts 
the acoustic wave passing from the source to each microphone. 

For the general 3-dimensional measurement geometry case, the set of equations in Amiet’^^Eq. A23 must 
be solved to locate the shear layer intersection point of a given ray. However, a 2-dimensional approximation 
can be used. Here, the 3-dimensional coordinate system is rotated about the flow axis such that the geometric 
line connecting the source and the observer lies within the new x — z plane. This rotation allows for the use 
of Eq. 1 from Amiet 21 under the assumption that the y intersection coordinate in the original coordinate 
system is co-planar with this rotated x — z plane. Alternatively phrased, the ray does not bend in any 
direction but the flow direction. The benefit is that the solution of a single equation nonlinear root-finding 
problem is significantly less expensive than the solution of a root-finding problem involving a system of non- 
linear equations. Additionally, this assumption automatically determines which shear layer is appropriate 
for each source-observer pairing by finding which shear layer intersects the new x — z plane. Testing of 
this assumption for various measurement configurations shows it to have minimal impact on the emission 
coordinates, ray direction and corrected levels with respect to the overall scale of the problem. 

This procedure outputs an acoustic level correction, as well as a model of the geometry of the acoustic 
ray path and emission coordinates. Amiet’s method corrects the data such that the emission path length is 
equal to the direct source-observer path length (as opposed to the refracted path length). Knowledge of the 
ray direction and emission path length is necessary for subsequent corrections, thus the shear layer correction 
is computed prior to the corrections for atmospheric attenuation and directivity. Emission coordinates are 
logged for use in the analyses reported in companion papers.®® 

IV. C. Atmospheric Attenuation 

For full-scale analysis of an HWB aircraft, model-scale data are desired up to 80 kHz. For source-observer 
distances on the scale of this test, atmospheric attenuation of the acoustic energy may be significant. A 
lossless acoustic “measurement” is often desired as an interim step for subsequent applications, such as 
model scaling and re-propagation calculations. A model of the attenuation loss may be applied to estimate 
this equivalent lossless propagated measurement. A standard for computing energy loss for a given frequency 
and distance as a function of atmospheric properties exist sp and is directly implemented and applied to 
the post-processed data. This standard provides the attenuation as a function of frequency, atmospheric 
pressure, temperature, and humidity. An example set of attenuation calculations is shown in Fig. [9] for a 
given atmospheric pressure and temperature. The plotted values have an uncertainty of ±10% at higher 
frequencies according to the standard.®^ The curves show attenuation in units of dB/inch, so at higher 
frequencies corrections of upwards of 10 dB may be applied to the data. The correction uses the emission 
path length provided by the shear layer correction calculation. 


IV. D. Microphone Directivity 


Microphones are not ideal signal observers, and have an influence on the recorded acoustic field determined 
by their installation characteristics.®^ The test was conducted with protective grids removed from the 
microphones (hence no grid correction is required), but the scattering microphone installation geometry 
must still be considered in the correction process. For free field microphone installation of the Briiel & Kjaer 
4138 1/8” pressure field microphones mounted along the sideline towers and overhead truss, a directivity 
correction table was provided by the manufacturer as a function of angle of incidence from the normal of 
the microphone face, 9, and frequency. This correction was applied with interpolation to the appropriate 
angles and frequencies. For most acquisitions, this correction was significantly greater than the corresponding 
pressure field correction for a given angle and frequency. 

For pressure field microphone installation, such as the Briiel & Kjaer 4938 1/4” pressure field microphones 
installed in the array plate, the majority of the installation influence is dealt with at the CSM construction 
phase when the signal level is divided by 2, approximating the installation as an infinite baffle. However, 
at higher frequencies, some wave averaging occurs over the face of the microphone diaphragm, and must be 
modeled. One simple way to model the microphone is as a flexible membrane installed in an infinite baffle. 
An analytic solution exists for the far held directivity of a radiating, simply supported circular membrane 
in an infinite baffle pi and is given by 


8 J 2 (kR sin 9) 
(kRsin9) 2 


( 14 ) 
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Figure 9. Example of atmospheric attenuation as a function of frequency and humidity, for atmospheric 
pressure of 101.3 kPa and temperature of 273 K. 


Here, k is the acoustic wavenumber computed using the speed of sound in humid airpSl ]{ is the membrane 
radius, and J 2 is the second order Bessel function of the first kind. Via reciprocity, the directivity pattern of a 
radiating, simply supported circular membrane should equate to the receptivity of a microphone diaphragm 


due to a source located at a given angle. Thus, the directivity function of Eq. ( 14 1 can be used to correct 
flush-mounted microphone levels by dividing the pressure value by D. It should be noted that as this 
formulation has a singularity at 9 = 0°, a refactoring may be preferable. Through a series of Bessel function 
identities Eq. (14 1 can be re-stated as 


4 1 

D (9) = J 0 (kR sin0) + - J 2 (kR sin 9) + - J 4 (kRsinO ) , 
o o 


(15) 


which, assuming access to good Bessel function libraries, avoids the computational trap for small 9. 9 is 
determined by the ray path computed in the shear layer correction calculation, in conjunction with the 
microphone orientation. As mentioned previously, for most acquisitions this correction was significantly 
smaller than the corresponding free field correction for a given angle and frequency. 

An example of the importance of these directivity corrections is shown in Fig. |10| This figure shows a 
comparison of autospectral measurements for the Broadband Engine Noise Simulators (BENS)P Here, the 
array plate was located such that the array center microphone was approximately normal to the axis of BENS 
exhaust, with no tunnel flow, and data were acquired for BENS operation. The acoustic hardware was then 
traversed such that the truss center microphone was located where the array center microphone had been 
located, and the measurement repeated. For the uncorrected data, recalling that the pressure doubling on 
the array face has already been removed, there is a significant difference between the two spectral estimates. 
However, when directivity effects are taken into account for both microphones, the spectra show agreement 
due to the significant shift of the truss microphone data from the free field correction. As expected based 
on the formulation, for a near-zero observation angle little correction occurs to the array microphone data. 


V. Summary 

An outline of the data processing methodology for the HWB acoustic test in NASA Langley’s 14- by 
22-Foot Subsonic Wind Tunnel is provided. The procedure for taking the data from the raw integer count 
stream recorded by the DAS to acoustic CSMs is shown as a step-by-step process, and discussed in terms 
of practical implementation. Existing code was updated from a legacy cluster Fortran code to a MATLAB 
program utilizing GPU computing, providing over an order of magnitude reduction in processing time and 
allowing for in-situ data analysis. Post-test modifications to the CSM construction process are reported. 
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Figure 10. Comparison of spectra for array and truss microphones approximately normal to the axis of 
BENS exhaust, without and with directivity correction. For this comparison, the free field correction 
applied to the truss microphone spectrum is the dominant factor in the data collapse. 


A novel transient rejection method is presented. The method analyzes a time series for transient bursts, 
and removes data blocks contaminated by the bursts. The method operates under the assumption that 
the desired signal to retain is a band-limited, stationary random process with lower block variance than 
the transient bursts to be removed. The method demonstrates the ability to separate contaminated blocks 
from uncontaminated ones. It is implemented as a completely automated process, with safeguards imposed 
should an aberrant signal show statistical behavior which causes the method to reject too much data. The 
method outputs CSMs, passed on to other analysis techniques, or autospectra, used with a baseline set of 
data corrections. 

Individual corrections to autospectral data are discussed, along with the assumptions involved in their 
use and the inherent limitations. Background subtraction can help remove facility background acoustics as 
well as electronic noise in the microphone measurements, under the assumption that the background noise is 
statistically independent from the acoustic sources of interest and stationary from acquisition to acquisition. 
Shear layer correction accounts for acoustic wave refraction through a shear layer by modeling the layer as a 
plane and solving an appropriate set of equations. A 2-D collapse of the problem reduces the set of equations 
to a single equation, accelerating analysis. Atmospheric attenuation reduces observed acoustic levels and is 
considered significant for the scales of this test. It is corrected by using an existing ANSI standard, under the 
assumption that the mean test section flow has little influence on the attenuation characteristics. Microphone 
measurement directivity corrections are performed using either supplied manufacturer data or an analytic 
model. A simple test case demonstrates expected collapse of the data between microphone installations when 
microphone directivity corrections are applied. 

The final product of these correction procedures is a set of lossless auto spectra corrected for observer 
interference, provided in emission coordinates at model scale. The data can then be used to perform a noise 
component analysis or the overall system-level noise assessment of the HWB N2A-EXTE. 
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